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What we (don't) know about cell 
mechanobiology and the nucleus 


Organized cell motion in an epithelial mono- 
layer (a single layer of cells) is important for 
tissue and organ function and homeostasis, 
but the mechanisms by which cells commu- 
nicate mechanical information to neighbour- 
ing cells in a coordinated and proportional 
way is poorly understood [1-3]. External and 
internal forces of the cell are transmitted 
through the cytoskeleton and can induce nu- 
clear deformation. Cell motility also requires 
dynamic, mechanical feedback between the 
cytoskeleton and nucleus for coordinating 
force transmission throughout the cell [4,5]. 
Yet, the mechanisms of nuclear positioning 
during cell migration are still not completely 
understood [6]. 


What we want to know 


In this study, we introduced a defective “stiff” 
cell within an epithelial cell monolayer with 
otherwise normal cells (Fig. 1 for details). 
This defective cell is the result of a mutation 
to the lamin A gene, which is an important 
structural element of the cell. The mutated 
gene is responsible for Hutchinson-Gilford 
Progeria Syndrome (HGPS), and contains a 
deletion of 50 amino acids, often called A5O 
lamin À (A5SOLA) or progerin [7]. Primary cells 
from patients with HGPS or cells exogenous- 
ly expressing ASOLA have stiffer nuclei [8], in 
addition to a host of other defects. Thus, we 
wanted to evaluate what effects emerge from 
the insertion of these “point defects” in a 
monolayer of otherwise normal cells. 
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Our Approach 


In this work, we utilized new temporal image 
processing methodologies to characterize 
nuclear morphology, motion, and shape dy- 
namics of normal nuclei in a monolayer with 
510% of NRK-52E cells expressing A5SOLA, 
These neighbouring nuclei (N) were com- 
pared with wild-type (WT) nuclei (monolay- 
er without the inclusion of ASOLA cells). We 
extended the analysis by investigating if any 
distance-dependency on the behaviour of N 
nuclei emerged. 


o Defected Nuclei 


€ Control Nuclei A " 
Region 7 VÁ .. etc .. 
“ Region 2 


O Region N II “a 
< Nuclear Centroid (E 
Cell , “a 


Figure 1- Subgroup division of N cells. Distance-dependency is 
evaluated by dividing the N nuclei into subgroups, according to 
the minimum distance to a AS0LA-expressing cell (red nuclei). 


Main Contributions and 
Discussion Highlights 


New methods from of nuclear 
tracking over time 


We showed that temporal fluctuations of nu- 
clear contours correlate with the stiffening 
of the nuclear lamina that is known to occur 
with the expression of ASOLA [8]. Other tem- 
poral techniques described in the full version 
of the paper can be used to determine cell 
mechanical properties (see Fig. 2 for com- 
putational analysis pipeline). Tracking the 
geometric nuclear centroid over time can 
approximate cellular motion accurately in 


monolayers of cells; finally, an intricate con- 
tour-based analysis of the nucleus allows for 
an insightful and detailed characterization of 
nuclear membrane fluctuations over time. 


(A) Contrast Normalization 
ImageJ 


Python 


Figure 2 - Overview of image processing methods used to ob- 
tain the intensity masks of each nucleus. A preprocessing step 
on Image) software includes contrast normalization of the 
image stack (A), which is then input to Icysoftware, where the 
remaining preprocessing takes place. Initial segmentation is 
achieved by the use of a Hierarchical K-means (HK-means) [9] 
method (B), which serves as reference for the application of an 
Active Contour plug-in [10], accomplishing accurate segmenta- 
tion and tracking of nuclei across the 21 frames (C). Finally the 
region outside of the nuclei masks is zeroed and an image stack 
containing only the nuclei of interest is obtained, which is used 
for processing and feature extraction (D). 


Stiffened inclusions influence whole 
monolayer structural properties 


Our findings revealed that cells are highly 
sensitive to the presence of mechanically 
impaired neighbours, leading to generalized 
loss of coordination in collective cell migra- 
tion (Fig. 3), but without seemingly affecting 
the potential for nuclear lamina fluctuations 
of neighbouring cells. Reduced translocation 
in neighbouring cells appears to be com- 
pensated by an increase in nuclear rotation 
and dynamic variation of shape, suggesting 
a “frustration” of cells and maintenance of 
motor activity. Interestingly, some character- 
istics of the behaviour of these cells appear 
to be dependent on the distance to a ASOLA 
cell, pointing to compensatory behaviour in 
response to force transmission imbalances 
in a monolayer. These insights may suggest 
the long-range impacts of single cell defects 
related to tissue dysfunction. 
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Figure 3 - Discoordination in cell motion is observed for the 
monolayer containing a ASOLA-expressing cell. A) Vector fields 
of nuclear centroids show direction and magnitude of cell dis- 
placement over time. The presence of a ASOLA-expressing cell 
leads to (B) overall lower cell speed, (C) increased changes in di- 
rection of trajectories and (D) lower migration efficiency which 
decreases nearly linearly with increasing distance to a mutated 
nucleus. Error bars represent 95% confidence intervals. 
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1. An unsolved problem 
of the 21st century 


Although rocket science, computer elec- 
tronics and biology have evolved to an out- 
standing level, the science in the protection 
of home buildings against wildland fires is 
Inexistent, even though the technology to 
promote its defence is available. This is even 
more unbelievable, since the cost of the re- 
search is very low and known knowledge in 
fire resistance is much more simple than 
rocket science or electronics. Therefore, tak- 
ing into to account the possibility to save 
thousands of lives every year, the cost-bene- 
fit analysis regarding the protection of hous- 
es for wildfires is positive. 


The main focus of this research project is to 
protect new dwellings from wildfires, using 
the concept of a fireproof envelope facade, 
with existent commercial non-combustible 
materials. 


2. From observation to an idea 


In 2017, two of the largest wildland fires 
(WLF) recorded in the history of Portugal 
(Pedrógão Grande and 15 de Outubro), had 
a major impact not only on the forest areas 
(burning 45.000ha+241000ha respectively) 
[Guerreiro et al. 2018; Viegas et al. 2017], but 
also on the communities, with a death toll 
of 114. Following these fires, new legislation 
was created [DL 76/2017 2017], only concern- 
ing measures related to the wildland-urban 
interface (WUI) (forest cleaning, vegetation 
height, distance from vegetation to dwell- 
ings, etc...). But interestingly, no extra mea- 
sure or care was taken, regarding the charac- 
teristics of the dwellings or the protection of 
the people caught by the fires. However, the 
analysis of expert opinions and questioners 
[Guerreiro et al. 2018; Tenreiro et al. 2018; 
Viegas et al. 2017], clearly showed that the 
dwellings that had burned, had constructive 
typologies without any passive anti-fire pro- 
tection. The disregard of this aspect, led to 
the burmned houses being rebuilt in the same 
way as the original, which means that they 
are likely to bum again in the next WLF. 
This important blank in protection mea- 
sures of new dwellings against WLF, led to 
the development of this project that aims to 
fill this gap in knowledge. To solve this prob- 
lem, fireproof construction guidelines (FCG) 
with incombustible materials are assembled 
and tested, for windows, doors, roofs and 
exterior walls. 


3. Wildland fire phenomena 


Global warming and subsequent climate 
change, has caused more often the occur- 
rence of wildland fires in the summer in the 
M editerranean, Europe, California and Aus- 
tralia, mainly in areas with large unattended 
forest areas and significant population occu- 
pation. In addition to the damage caused in 
the forest, these fires cause the destruction of 
dwellings and occasionally lead to the death 
of residents, therefore, it is necessary to im- 
plement constructive preventive measures. 


a) Masonry Dwellings 


Genet A 


c) Wood Roof Structure 


228 corpo > , E 
d) Masonry Dwellings with Wood Roof 
Figure 1 - Example houses that burned, adapted from Portu- 
guese press SIC, RTP and DN. 
The narrative of the destructive power of 
fire on nature itself was reinforced during 
the 20th century. The intensity and hazard 
of WLF, depends on 3 mechanisms of fire 
spread, which are fuel, wind and topography 
[M ueller 2016]. M oder society has generally 
taken two approaches, forest management 
[Davis 1979] and suppression [Pyne et al. 
1996] to dealing with the unavoidable exis- 
tence of wildland fire. 


Despite the fact that both management 
and suppression tactics have existed for a 
long time, the ideal means of implement- 
ing these strategies, particularly in terms of 
management, has yet to be achieved [Ste- 
phens et al. 2005). 


As human activity and development has 
pushed ever further from urban centers, the 
past few decades have seen the growth of 
what is often referred to as the Wildland Ur- 
ban Interface (WUI) [Glickman et al. 2001]. 
The growth of the WUI becomes problematic 
when we consider that it can be directly as- 
sociated with settlements and infrastructure 
in areas which are particularly risk prone for 
WLF [Radeloff et al. 2005]. 


There are a number of tools available to help 
address and manage the WLF problem, as 
defined above. At the moment very small rec- 
ommendations exist concerning construction 
guidelines, to ensure that any new dwelling 
possesses enough defences against a WLF. 
This project plans to fillthis gap in knowledge. 


4. Main objectives 


The main objective of this project, is to study 
the efficiency of the FCG with incombusti- 
ble materials (existing on the market), after 
a WLF as occurred, regarding that the dwell- 
ings don't suffer any damage, when hit by a 
WLF. To analyse the efficiency of FCG in post 
fire conditions, regarding its fire resistance 
and reaction, these will be subjected to ther- 
mal standard fire tests [Razdolsky 2022]. 


The main objectives are: 


1) Perform a fire experimental campaign for 
specimens and reduce models, for incom- 
bustible materials and fireproof construc- 
tive solutions, using standard and wildfire 
fire curves, in order to check the level of fire 
resistance and reaction of the initial chosen 
solutions. 


2) Execute parametric numerical fire cam- 
paign, in order to extrapolate the results of 
goal 1) to other structural or geometric ar- 
rangements, maintaining the level of fire resis- 
tance of the previous experimental campaign. 


3) The Proposal of a construction manual, to 
be distributed by the municipalities, in resi- 
dential areas with a high risk of WLF, inform- 
ing the population about the construction 
fireproof solutions to be adopted to protect 
their homes from an external fire, depending 
on their exposure conditions. 


5. Expected results and 
contribuition to scentific and 
technological knowledge 


The work developed, will be a very import- 
ant innovation regarding passive protection 
of dwellings from WLF in Portugal, allow- 
ing the release of the firefighter to fight the 
fire directly in the forest field. This research 
project will promote the increase of safety in 
dwellings, population and human resources 
used in WLF. 


Beyond that, it will allow to reduce the cost of 
repairing or rebuilding dwellings after WLF. 


It is then possible to; 


1) The reduction of human drama due to 
the scourge of fires. Besides saving lives, as 
added value this solution prevents a home 
building from being consumed, giving the 
individual the possibility of not losing what 
it took years to build. In this case, the affec- 
tion for the house, memories and everything 
that is inside. 


2) À cost reduction during the post-fire, since 
it wouldn't be necessary any rehabilitation 
or reconstruction of dwellings. When com- 
pared with the fire in Pedrogão Grande in 
2017, 54 residential houses were reconstruct 
and 102 were rehabilitated, with a total cost 
of 9.9 million euros [Pinto 2018]. It is esti- 
mated, that the average cost of construction 
with FCG is around 250€/m2, and the cost 
of rebuilding a dwelling is above 1000€/ m2. 


3) Firefighters have more flexibility in resi- 
dential areas, allowing the firefighters to ef- 
fectively contain and control/extinguish the 
fire more effectively in the forest field. 


6. Initial results 


The design of an “fireproof house” implies 
that it: (1) must withstand the near passage 
of a WLF; (li) it must be protected against the 
entrance of firebrands; (lil) and in urban situ- 
ations, should be protected against the pos- 
sible fire of nearby dwellings. At the moment, 
with less than 1 year of research, the team as 
assemble and published fireproof construc- 
tive guidelines [Arruda et al. 2021], using ma- 
terials that exist in the construction industry, 
and that are also common known for their ex- 
cellent thermal behaviour. 


a) Fireproof roof solution with glass or rock wool 


c) Fireproof wall ETIC with rock wool 


Figure 2 - 
and walls. 


Initial fireproof solutions for roofs, window, doors 


Contrary to normal ISO fire curves, a WLF 
possesses a maximum temperature of 800 
to 900ºC for 3to 10 minutes period, depend- 
ing on wind, fuel and topography [Pyne et al. 
1996]. In addition, the maximum tempera- 
ture of a firebrand is above 400ºC for more 
than an hour (Thomas, 2018 4207). There- 
fore, the ISO fire curves are more severe than 
a natural wildfire. For this reason, new stan- 
dard wildfire curves were proposed [Cantor et 
al. 2022], in order to numerically simulate the 
action of a WLF in the dwellings exterior, and 
therefore assembling a more economically 
fireproof solution. 


Partners were also found in the construction 
industry to collaborate in this research proj- 
ect, in order to promote the test of new eco- 
logic with a negative carbon footprint, that 
could be used as fireproof solutions. 


7. Project Web Page 


More information on the outputs, and deliv- 
erables of this research project can be found 
in the following link: whp.tecnico.ulisboa.pt 
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Palavras-chave: Internet das Coisas, Lumi- 
nescência, Fotónica, Funções Físicas Não 
Clonáveis 


A rápida expansão da Internet das Coisas in- 
crementou a necessidade de métodos rápi- 
dos efiáveis para a identificação de objetos e 
processos. Neste trabalho, propomos iden- 
tificadores fotónicos, compactos, baseado 
em funções físicas não clonáveis, e adequa- 
dos à utilização com recurso a dispositivos 
móveis (smartphones). Estes identificadores 
físicos, de baixo custo, são produzidos com 
materiais híbridos orgânicos-inorgânicos lu- 
minescentes flexíveis e dopados com iões 
lantanídeos, apresentando pares desafio-re- 
sposta únicos e compatíveis com procedi- 
mentos de autenticação de dois fatores. 


Este trabalho poderá ser consultado na ín- 
tegra na seguinte publicação [1]. 


Introdução 


A Internet das Coisas (conhecida pelo acrón- 
imo loT, da designação em língua inglesa 
“Internet of Things”) é uma rede de objetos 
conectados, interagindo entre si em tempo 
real e sem supervisão humana, sendo fun- 
damental no desenvolvimento das futuras 
cidades inteligentes [2]. A disseminação em 
larga escala desta rede exige a segurança dos 
dispositivos conectados através de soluções 
de autenticação que garantam e validem a 
identidade dos objetos interligados. De en- 
tre as soluções de segurança propostas para 
dispositivos loT, destacam-se as funções físi- 
cas não clonáveis (conhecidas pelo acrónimo 
PUFs, da designação em língua inglesa de 
“Physical Unclonable Functions”), por con- 
jugarem singularidade e irreprodutibilidade. 
As PUFs são soluções inovadoras e de baixo 
custo, baseadas em perturbações físicas ale- 
atórias, caracterizadas por apresentarem 
uma resposta complexa a um determinado 
desafio a qual é transformada numa “im- 
pressão digital” única [3] 


O conceito de PUF foi introduzido em 2001, 
com recurso acircuitos eletrónicos [4] Porém, 
essa implementação no domínio eletrónico é 
vulnerável a ciberataques, especialmente aos 
promovidos com métodos de aprendizagem 
automática [5]. Para colmatar essa vulnerabi- 
lidade, surgiram PUFs no domínio ótico que 
geram respostas mais complexas e que ofere- 
cem uma maior resiliência à manipulação [6]. 
O princípio de funcionamento subjacente às 
PUF óticas baseia-se numa fonte Ótica que 
ilumina (interroga) um identificador físico, 
sendo o sinal refletido, espalhado, dispersa- 
do ou transmitido, a resposta que contém 
a informação unívoca da identificação [7,8]. 


Considerando que os smartphones são uma 
tecnologia omnipresente, que incorpora dis- 
positivos de captação de imagem capazes de 
ler as respostas óticas de PUFs, com capaci- 
dade computacional capaz de as processar, 
e com uma interface de comunicação ade- 
quada à ligação a uma rede externa, podem 
contribuir significativamente para a leitura da 
informação contida nas PUFs [9]. 


Neste trabalho, reportamos a produção de 
PUFs fotónicas funcionais, flexíveis, produz- 
idas a baixo custo, com métodos escaláveis 
e energeticamente eficientes, utilizando um 
híbrido orgânico-inorgânico dopado com 
iões de lantanídeos (Ln*t=Eu+Fig.1). As 
PUFs são materializadas através de identifi- 
cadores físicos que apresentam uma miríade 
de respostas quando interrogados por um 
sinal Ótico gerando pares de desafio-respos- 
ta (conhecidos pelo acrónimo CRP da des- 
ignação em língua inglesa de “challenge re- 
sponse pair”). Propomos 2 tipos de CRP, 1) 
baseados na interação com sinais Óticos co- 
erentes, gerando padrões de granitado laser 
(speckle), (Fig.Lb), que podem ser adquiri 
dos pela câmara de um smartphone; e ii) ba- 
seadas em sinais Óticos incoerentes que pro- 
duzem padrões de luminescência, analisados 
por microscopia hiperespectral (Fig.1c). Esta 
dualidade aumenta a complexidade da res- 
posta e adiciona uma camada adicional de 
proteção, permitindo um esquema de auten- 
ticação de dois fatores (granitado e emissão). 


Reportamos uma prova de conceito, conce- 
bida para um cenário real, onde os identifi- 
cadores físicos desenvolvidos são aplicados 
a um cartão bancário, demostrando-se que 
estes identificadores possuem as característi- 
cas mecânicas e óticas adequadas à imple- 
mentação em tais cenários, sem comprome- 
ter a metodologia de autenticação [1]. 


Metodologia 
Produção e caracterização dos identificadores 


Foram produzidos identificadores físicos, 
baseados no precursor orgânico-inorgânico 
híbrido não hidrolisado, d-UPTES(600), e o 
complexo Eu(tta)32H 20 (Hita=2-tenoyltri- 
fluoracetona), sendo o material resultante 
denominado como dU6Eu [1]. Os híbridos 
podem ser produzidos numa variedade de 
formas, adaptadas a diferentes aplicações, 
como ilustrado no caso do cartão bancário 
(Fig.Le). Os defeitos e irregularidades obser- 
vados na imagem microscópica Ótica (Fig.Ld) 
mostram a microestrutura responsável pela 
aleatoriedade das respostas. 


Para a validação da metodologia aqui apre- 
sentada, foram adquiridas imagens de granit- 
ado laser de dois identificadores físicos difer- 
entes (CRPLA e CRP LB). Para os mesmos 
identificadores físicos, as CRP baseadas nos 
padrões de emissão luminescentes adquiri- 
das estão reportados em [1]. 


Avaliação da autenticação 


A extração da chave de identificação a par- 
tir das imagens é realizada através de uma 
função de hashing, baseada na transformada 
discreta de cosseno (conhecida pelo acróni- 
mo DCT da designação em língua inglesa de 
“discrete cosine transform”), nomeadamente 
a DCT 2D tipo Il, equação (1) [10]. 


C; (u, v) = 


(o + 2) e + 2) (1) 
os[————— |cos|————., 


2m 2m 


As imagens são divididas em mxm pixels, 
não sobrepostos, e a DCT é aplicada a cada 
pixel B;(),k), resultando no coeficiente C(u,v) 
[10]. Os termos superiores esquerdos da 
matriz DCT, utilizados neste trabalho, corre- 
spondem aos termos de baixa frequência na 
resposta baseada na imagem, enquanto os 
coeficientes inferiores direitos estão associa- 
dos aos termos de alta frequência [11]. 


Após a obtenção das chaves de identificação 
da imagem a autenticar, é possível correlacio- 
nar com os valores previamente armazena- 
dos eos resultados são utilizados para avaliar 
o desempenho da PUF e determinar as carac- 
terísticas do processo de autenticação. A fig- 
ura de mérito mais utilizada para quantificar 
o desempenho das PUF, e adotada neste tra- 
balho, é a distância de Hamming normaliza- 
da (HD) [2]. As distâncias intra-Hamming 
(intra-HD) e inter- Hamming (inter-HD) são 
calculadas para comparar as chaves extraí- 
das, respetivamente, do mesmo identificador 
e de identificadores diferentes [13]. Foram ex- 
traídos dos identificadores chaves de difer- 
entes tamanhos e comparados para os aval- 
lar assim como à sua utilização em termos 
de singularidade, fiabilidade, uniformidade, 
bit-aliasing [14], e probabilidade de erro na 
identificação [15]. 


Resultados e discussão 


O espectro de absorção do dU6EU (Fig.1f) 
revela uma ampla faixa de absorção na região 
espectral UV/azul, assegurando a transpar- 
ência sob a luz solar. Sob radiação UV, exibe 
uma cor vermelha das transições intra-4f6 
(Fig.Lg). A capacidade de controlar a cor e 
transparência dos identificadores através do 
ajuste do comprimento de onda da radiação 
incidente prevê a possibilidade de definir 
diferentes CRP (CRP, Fig.1Lb) e emissão 
(CRP2, Fig.Lc). Estas duas regiões espectrais 
representam a gama de operação do PUF, no 
visível (480-800 nm) eUWV/ azul (250-450 nm), 
respetivamente, para o granitado e emissão. 
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Figura 1 a) Esquema da autenticação de dois fatores utilizan- 
do PUFs fotónicas. O identificador exibe pares desafio-resposta 
distintos que permitem a autenticação de dois fatores basea- 
da em b) granitado laser quando irradiado com luz coerente 
e c) emissão sob radiação UV (LEDs), suportados pela micro- 
estrutura do identificador baseada em defeitos observados na 
imagem de d) microscopia ótica; espectros de f) absorção, 9) 
excitação (azul) e emissão (vermelho), e) fotografias do identi- 
ficador físico embutido num cartão bancário sob iluminação [1]. 


Sob o mesmo desafio, foram consideradas 
as respostas de granitado de dois identifi- 
cadores, denominados CRPLA (Fig.2.a) e 
CRPLB (Fig.2.b). A avaliação dos identifica- 
dores está representada nas Fig.2.c e d. A 
comparação entre os identificadores de 28 
bits através do HD é mapeada na Fig.2.e, 
onde é possível localizar as zonas azuis que 
correspondem a baixos valores HD (perto de 
0) entre as respostas do mesmo identificador 
(intra-HDs paraCRPLA eCRP LB), easzonas 
vermelhas com valores H D superiores (perto 
de 0,5) que correspondem à comparação en- 
tre as respostas dos dois identificadores dis- 
tintos (inter-H Ds). Os histogramas para um 
elevado conjunto de respostas podem ser 
descritos por duas distribuições Gaussianas 
obtidos das H Ds ajustadas a funções de den- 
sidade de probabilidade (Fig.2.f), idealmente 
centradas a O (intra-HD) e 0,5 (inter-HD). 
Foi possível obter um limiar de decisão óti- 
mo (dopt) que minimizou a probabilidade 
de erro de identificação, com donj= 0,31 
para um probabilidade de erro de aproxima- 
damente 10-!!, inferior a valores previamente 
reportados [16]. 
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Figura 2. Imagens de padrões de granitado laser a) CRP LA eb) CRP LB. c),d) Parâmetros de avaliação da PUF. e) Mapeamento das 
HD calculado a partir das respostas de 28 bits, onde as regiões azuis correspondem à comparação entre as respostas adquiridas do 
mesmo sinal (CRP LA vs CRP LA, números de imagem 1a 50 e CRP ÉB vs CRP LB, números de imagem 5la 100) e as regiões vermel- 
has correspondem à comparação entre CRP LA vs CRP É:B. f) Histograma dos valores HD de 100 respostas com base numa chave 
de 28 bits. As funções de densidade de probabilidade intra-H D e inter-HD têm um valor médio de 0,06 e 0,50 e, respetivamente, 


um desvio padrão de 0,04 e 0,02 [1]. 


A viabilidade de adicionar um segundo par 
de autenticação independente (CRP2) para 
cada ficha baseia-se na fotoluminescência 
dos dispositivos, caracterizada por um pa- 
drão baseado na intensidade aleatória das 
emissões. Os detalhes desses resultados 
podem ser consultados na referência [1], 
mostrando também uma implementação 
promissora utilizando as características de 
luminescência. 


Conclusão 


A possibilidade de utilizar um smartphone 
para autenticação permitirá traçar o caminho 
da integração de PUFs óticas com dispositi- 
vos eletrónicos, abordando um dos aspetos 
cruciais no campo da autenticação ótica. 
Foi demonstrada a utilização de dispositivo 
de autenticação móvel funcional utilizan- 
do PUFs fotónicas baseadas em materiais 


híbridos orgânicos-inorgânicos selecionados 
dopados com iões Eu?* (dUGEU). O mapea- 
mento HD e os resultados dos histogramas 
apresentaram informações valiosas sobre os 
parâmetros de autenticação a serem utiliza- 
dos, nomeadamente, o limiar de decisão óti- 
mo e a probabilidade de erro para o sucesso 
da autenticação. A caracterização das chaves 
de autenticação extraídas das respostas mos- 
trou a elevada uniformidade (50%), fiabi- 
lidade (>90%), unicidade (>40%) e a baixa 
probabilidade de erro (<10-1). Este sistema 
oferece um paradigma completamente novo 
na definição de chaves de identificação at- 
ravés de imagens adquiridas por um smart- 
phone que podem ter diversas áreas, uma 
vez que proporciona um processo de autenti- 
cação remota fácil, barato e rápido. 
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1 Introduction 


Throughout history there are different 
examples of different organizational 
structures, ranging from strictly hierarchical 
organization such as that embodied by 
the military to fully-distributed consensus 
based decision-making processes used 
in some worker cooperatives, and many 
different structures in between. Despite 
this wide variety of existing govemnance 
structures, current cybersecurity tools focus 
on assuming governance models that are 
strictly hierarchical, focusing on singular 
administrators with large amounts of power 
over the system. 


Though this model of system administration 
works well for hierarchical organizations 
such as military, traditional companies, and 
others, it does not address the needs of all 
organizations. More democratic entities 
such as activist groups, worker cooperatives, 
some not-for-profit organizations, trade 
unions, and more tend to run themselves 
more democratically, avoiding or attempting 
to avoid the centralization of power into 
the hands of a few actors. However, current 
cybersecurity tools and techniques often 
assume a small set of administrators, causing 
such a centralization of power to occur. 


In my work, | have been reimagining 
cybersecurity from the ground up to include 
tools and techniques that better fit more 
horizontal organizations. This includes 


re-imagining access control to allow for 
collectively-managed digital assets, using 
systems built on checks and balances to 
protect against rogue administrators, and 
re-imagining threat modeling protocols 
to consider how a given tool or technique 
could cause power centralization within an 
organization. 


Though this long-term goal may seem 
ambitious or overwhelming, work has already 
started in this direction. In the following two 
sections | describe the work that has already 
been done in the area, as well as the currently 
ongoing work. In Section 4 | discuss the 
future vision of this work, both in projects 
that will begin shortly as well as longer- 
term visions of how these solutions will fit 
together to hopefully cause a paradigm shift 
in cybersecurity. 


2 A First Step in 
Horizontal Security 


During the Occupy M ovement, some admin- 
istrators of non-hierarchical activist groups 
began locking their fellow members out of the 
group's social media pages and posted sever- 
al harmful posts about the activist group they 
belonged to [5]. These administrators used 
their privileged access to the group's pass- 
words to centralize the group's communica- 
tion capabilities, bypassing collective control. 
Although it may be tempting to think of this 
as a one-off problem experienced by a single 
activist group, Gerbaudo demonstrates that 
this problem is an emerging issue among 
many activist groups [3]. 


Upon first glance, the obvious answer to this 
seems to be some secret sharing scheme 
- Sharing the password among multiple 
entities so that no one user has the whole 
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password, and posts can be made only with 
permission of multiple users. However, fur- 
ther inspection shows this solution to be in- 
sufficient; a rogue administrator could sim- 
ply use their administrator authorityto cause 
the centralization of power elsewhere in the 
organization's technological infrastructure, 
such as email servers or communication 
channels such as Signal or Telegram groups. 
As such, we decided to revisit the adminis- 
tration problem itself, specifically as it ap- 
plied to access control. 


We attempted to resolve this issue by creating 
a form of access control focused on the collec- 
tive and its resources called COLBAC [2]. This 
system focuses on ensuring that all access 
to a resource owned by a collective such as 
an activist group, trade union, or some other 
democratic group is granted by a democratic 
process, and that no action taken by a mem- 
ber of the collective could cut the collective 
off from their resources. This solution worked 
by using a voting system that allowed users 
to individually vote on a given action, like a 
referendum in modern governance schemes, 
or to vote on electing an administrator to per- 
form limited set of actions, like an election in 
modem representative governance schemes. 


Unlike modem governance schemes, howew- 
er, any elected administrator can be removed 
from power at any time with a successful ref- 
erendum. Thus, if an administrator abuses 
their power to attemptto disrupt the goals of 
the collective, the collective can regain their 
access and revoke access to the administrator 
quickly and effectively. This is implemented a 
logical division of a system or network into 
three or more spheres, the collective sphere, 
which holds the collective property of the or- 
ganization and all of the system critical files 
and data, the user sphere(s), which hold files 
that belong to individual users and are not 
the property of the collective organization, 
and the immutable sphere, which holds logs 
necessary to ensure transparency on behalf of 
the administrators. Any action to modify a file 


in the collective sphere would require a token 
issued bythe collective by some collective de- 
cision-making process. Due to space limita- 
tions we omit many details here, but interest- 
ed readers can refer to the original paper [2] 
for more details. 


Unlike the previously considered solution of 
secret sharing, this solution is much more 
general, and applies equally to administering 
email servers as it does to posting on social 
media. Like all systems, however, COLBAC 
is imperfect. It is important to note that this 
solution does not prevent the issue from oc- 
curring in the first place, but rather allows a 
democratic organization to be resilientto this 
form of attack and recover access to their re- 
sources quickly after this form of attack has 
taken place. More, there are suspected us- 
ability issues with the system design, and the 
general approach means that while it may ap- 
ply well to many different situations, it may 
be too general to fit any of those situations 
perfectly. As such we decided to move for- 
ward with future work that took the spirit of 
COLBAC and applied to specific situations in 
order to create a system that fits the needs of 
a given specific application. 


3 Current Work in 
Horizontal Security 


In order to further explore the concept of hor- 
izontal security and the various ways in which 
it can be achieved and applied, we decided to 
move our exploration into multiple different 
directions. These directions included study- 
ing DAOs, or decentralized autonomous or- 
ganizations and the technologies they used, 
as well as working on threat modeling for 
organizations with less-hierarchical organiza- 
tional structures, and on creating an access 
control system for Parliament administrators. 
Due to space limitations, here we describe 
only the ongoing work in the latter two cases. 


3.1 Considering Governance 
Structure in Threat Modeling 


In our COLBAC work we made the assump- 
tion that the hierarchy assumed by certain 
cybersecurity tools and techniques, such as 
administration techniques and access con- 
trol techniques, impose hierarchy on the 
organizations that use it. Though there ap- 
pears to be plenty of anecdotal evidence to 
suggest this being correct [8, 5, 3, 4], to our 
knowledge there has been no definitive way 
of measuring the influence the use of a tool 
or technique has on the governance structure 
of an organization. Therefore, in this work we 
are proposing a threat modeling protocol for 
organizations to determine the threat of cen- 
tralization of their governance scheme based 
on the techniques and technologies they use, 
similar to the way that current threat mod- 
eling tools like STRIDE [1] help developers 
determine what attacks their software could 
fall victim to. This threat model protocol 
will examine the governmental structure of 
an organization and the technologies and 
techniques they use, and will output a set of 
threats that could force changes in the gov- 
ermance system based on abuse of the tech- 
nologies or techniques they rely on. This is 
ongoing work in its early stages in joint col- 
laboration with Sebastián Chimal M ontes de 
Oca from IST, Santiago Torres-Arias, an as- 
sistant professor at Purdue University, and 
Jessica Feldman, an assistant professor at 
American University of Paris. 


3.2 Defendin 
from Insider 


Parliament 
track 


Perhaps the most obvious and iconic exam- 
ple of an organization that is meant to be 
non-hierarchical is the democratic institution 
of Parliament. Though some parliamentari- 
ans do have different responsibilities of vary- 
ing importance, when it comes to voting on 
legislation each parliamentarian is meant to 
have an equal vote and equal power. More, 
Parliament is also a very high-risk organiza- 


tion. Itis made up of many different political 
parties that sometimes cooperate and some- 
times compete in order to pass the political 
agenda most closely aligned to their party's 
political ideology. Though sometimes these 
ideologies align, often times these parties” 
ideologies conflict, and lead to more compe- 
tition than cooperation. 


Currently, Portuguese parliamentarians rely 
on a shared infrastructure hosted by the Por- 
tuguese parliament itself. Inside Parliament 
there is a team of trusted administrators who 
oversee the shared resources, which include 
networking equipment, computer systems, 
email services, and more. Currently these 
administrators have large amounts of power 
over the network - after all, they need to fix 
any issues that may arise in the infrastructure 
- and have very little oversight into how they 
choose to use it. Consider, for example, that 
an administrator in Parliament that is respon- 
sible for the email services may have a very 
strong political leaning and support a specific 
political party. Given that this administrator 
has nearly unrestricted access to the email 
system, this administrator may be able to use 
this access to spy on the parties that compete 
with their favored party. Worse still, this at- 
tack could be silent, as the rogue administra- 
tor could also remove or modify any logs that 
demonstrate that this access was performed. 
We seek to prevent some attacks that could 
be performed by rogue administrators as well 
as provide resiliency to other attacks by creat- 
ing a COLBAC-like system to provide checks 
and balances within parliament infrastruc- 
ture administration. Instead of allowing ev- 
ery member of the system to vote, such as in 
COLBAC, in this instance we propose an ad- 
ministrative council that performs COLBAC 
operations. This council is made up an equal 
number of representatives from each political 
party ensuring that any action against that 
political party by a proposed administrative 
action is known to that party and that the par- 
ty has a say against it. This council will then 
go through all of the steps of COLBAC. 
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A prototype of this system is currently under 
development. The prototype uses a Byzantine 
Fault Tolerant Consensus [6], modified PAM 
modules, and asymmetric cryptography to 
create, distribute, and validate COLBAC to- 
kens. After implementation, we will validate 
the prototypes effectiveness in terms of ef- 
ficiency and security described as both pre- 
venting some types of rogue administrator 
misconduct, such as performing unautho- 
rized actions on the system, and providing 
resiliency to other forms of rogue administra- 
tor misconduct, such as abusing legitimate 
access to attack parliament or the political 
parties that make it up. 


Focusing on parliament administrators al- 
lows us to focus on a practical implementa- 
tion of the ideas of COLBAC with a specific 
goal in mind, which helped in solving some 
of the more practical questions of using 
COLBAC in a distributed system or network 
setting. For example, in the general version 
of COLBAC discussed in [2] we did not fo- 
cus on the details about how distributed or 
complex systems could use COLBAC by rep- 
licating tokens and logs across the an entire 
organization's Infrastructure. By focusing 
on a prototype implementation of a COL- 
BAC-like system in a complex infrastructure 
like Parliament's computer and networking 
systems, we were able to build a system that 
addressed that concern efficiently using BFT 
consensus. 


This is joint work with Margarida Espanhol 
Lopes and Professor Miguel Correia here at 


Instituto Superior Técnico, and is nearing 
Its final stages. We will enter in contact with 
Parliament soon to attempt to gather their 
reaction to our prototype work and use 
their feedback to improve our solution to fit 
their needs. 


4 Future Work 


Similar to our work on developing a more 
horizontal security tool for parliament, future 
work will see us applying the idea of horizon- 
tal security to organizations and groups that 
rely on democratic governance structures. 
Towards these ends we are beginning discus- 
sions with horizontally run activists groups, 
free software movements, and more in vari- 
ous countries, including M éxico, the United 
States of America, Portugal, and more. We 
will use these conversations to inform the 
creation of more tools for these organiza- 
tions, and develop individual prototypes and 
software to fit their security needs. Examples 
of these needed tools could include commu- 
nity moderation for instant messaging, chat, 
and social media services, democratically in- 
cluding non-human agents such as animals 
or plants into technological decision that af- 
fect them, and more. When these individual 
projects are done we will take what we have 
learned from all of them and create a suite of 
horizontal security tools meant to fit a wide 
array of democratically run organizations, po- 
tentially ending with a new operating system 
on which they can host their websites, email 
servers, and more. 


5 Conclusions 


In his work The Moral Character of Cryp- 
tographic Work, Rogaway discusses the ten- 
dency of cryptographic research (and, | argue, 
general security research) to benefit organi- 
zations and people that already have power 
[7]. A part of this is the focus that security 
tools make on fitting the needs of hierarchi- 
cal organizations such as large businesses 
and military organizations while neglecting 
to research tools that focus on less powerful 
communities such as citizen assemblies, ac- 
tivist groups, worker, consumer, and housing 
cooperatives, trade unions, and more, which 
organize themselves in a more democratic 
manner. However, in our recent and future 
work we are working on expanding the vi- 
sion of the cybersecurity research communi- 
ty to include these groups. In this article we 
discussed our first steps into this world: the 
development of COLBAC, and our current 
work focusing on threat modeling for gov- 
ernance structures and implementation of a 
COLBAC-like system for Portuguese parlia- 
ment. In future work we will be working more 
closely with activist groups, free software 
movements, and more to develop more se- 
curity tools and techniques specifically build 
for their needs and organizational structures. 
By doing so, we can help these organizations 
remain secure, perhaps aiding their growth 
until they can rival the size of their business 
and military counterparts. 
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Abstract 


Cataract, highly prevalent among old popula- 
tion worldwide, is generally treated by replac- 
ing the opacified crystalline lens by an intra- 
ocular lens (IOL). Although this procedure 
is safe and effective, there are some com- 
plications that can occur after surgery, like 


posterior capsule opacification (PCO) and 
endophthalmitis, sometimes with serious 
consequences. IOLs with different designs 
and materials have been studied to reduce 
PCO. For the prevention of endophthalmitis, 
intracameral injection followed or not by topi- 
cal treatment with antibiotics and anti-inflam- 
matories are usually prescribed. The objective 
of this work was to investigate the possibility 
of using IOLs as controlled release platforms 
of two drugs, the antibiotic moxifloxacin 
(MXF) and the anti-inflammatory ketorolac 
(KTL), to advantageously substitute the usu- 
al treatment. Two types of IOLs were chosen, 
hydrophobic and hydrophilic. Hydrophobic 
IOLs have shown better results in the pre- 
vention of PCO because they adhere better to 
the posterior capsular bag, while hydrophilic 
IOLs are advised in the case of patients with 
uveitis, glaucoma or diabetes. The IOLs were 
loaded with MXF+KTL and sterilized by high 
hydrostatic pressure. Both IOLs reduced the 
tendency for adhesion of lens epithelial cells 
comparatively to unloaded IO Ls. In vivo tests 
were done to compare the drug's concentra- 
tion in the aqueous humor after eye drops ad- 
ministration and drug-loaded IOLs implanta- 
tion. The developed IOLs released MXF and 
KTL at therapeutic levels, in a sustained way, 
which contrasts with the eye drops prophylax- 
is. Neither ocular infection/ endophthalmitis 
nor PCO signs were detected. This article is a 
short version of our publication J. Controlled 
Release 326 (2020) 245-255, where all exper- 
imental details and full experimental results 
can be found. 
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Introduction 


Cataract is the leading cause of blindness 
worldwide, affecting —94 million people [1]. 
The huge medical, social and economic im- 
pact of this disease justifies the increasing 
number of studies on the prevention and 
treatment of related complications. In fact, 
although cataract removal is a common- 
ly performed surgery, it is not free of risks. 
The posterior capsule opacification (PCO), 
which may occur months or years after sur- 
gery decreasing significantly visual acuity, 
has an incidence that can reach 50% [2]. In 
contrast, post-cataract surgery endophthal- 
mitis is a rare (0.04-0.2%) but severe com- 
plication, and usually leads to limited visual 
recovery [3,4]. 


PCO results from the adhesion of lens epithe- 
lial cells (LECSs) to the surface of the IOL fac- 
ing the posterior capsular bag. The reduction 
of the PCO can be addressed in three ways: 
adequate choice of the IOL biomaterial [5,6], 
design of the lenses [5,7,8] and drug vehicu- 
lation using the IOLs as delivery devices [9- 
11]. Hydrophobic polyacrylic IOLs have been 
recognized as being associated with a lower 
PCO incidence than hydrophilic polymethyl- 
methacrylate and silicone IOLs [5]. They can 
adhere more easilyto the collagen membrane 
of the posterior capsular bag and minimize 
the possibility of LECs reaching the lens sur- 
face [2]. However, this type of lenses has a 
high tendency to adsorb extracellular matrix 
proteins and inflammatory cells, which may 
lead to other complications, such as iris pos- 
terior synechiae and anterior capsular opaci- 
fication. For this reason, although hydro- 
phobic IOLs are the most commontly used, 
they are not recommended for patients with 
uveitis, glaucoma and diabetes who receive 
hydrophilic IOLs instead [13] Other authors 
claim that PCO incidence is independent of 
the biomaterial's nature [14], thus the role of 
the biomaterial remains a controversial is- 
sue. Concerning the lenses design, reduced 
PCO incidence was found in IOLs with sharp 


edges [5-7]. Finally, in recent years, the use of 
IOLs as drug delivery devices has been iden- 
tified as a potential strategy for PCO prophy- 
laxis [9-11,15]. 


Endophthalmitis is mainly caused by 
gram-positive bacteria such as Staphylococ- 
cus aureus, Staphylococcus epidermidis and 
Streptococcus pneumonia and may assume 
an acute form that occurs in the first 48 h fol- 
lowing surgery or appear more than 6 weeks 
afterwards. Furthermore, cataract surgery 
may induce an inflammatory response, lead- 
ing to discomfort, compromised vision and 
cystoid macular edema (CME) [16]. To pre- 
vent or treat these complications, eye drops 
of antibiotics and anti-inflammatories are pre- 
scribed to be administered during a relatively 
long period of time and with a very frequent 
posology [17,18]. Topical instillation leads to a 
low drug bioavailability (<5%), dueto the rap- 
id clearance of drug from the eye surface [19], 
and to adverse side effects due to absorption 
at the systemic level [20]. Intracameral injec- 
tions, which are also an option for antibiotics 
[21], are rarely used, despite the observed re- 
duction in the endophthalmitis rate [22] due 
to toxicity concerns. Recent studies point out 
to the possibility of using IOLs as the vehicle 
to deliver, in a sustained way, antibiotics and 
anti-inflammatories for endophthalmitis pro- 
phylaxis [23-26]. 


In the present work, the two main problems 
associated with IOL implantation (endoph- 
thalmitis and PCO) were tackled and possi- 
ble solutions were tested through in vitro and 
in vivo tests. Acryic based hydrophobic (GF) 
and hydrophilic (CI26Y) lenses were used. 
Hydrophobic lenses, which have a low water 
content, were modified through the coating 
of the anterior side with a hydrophilic layer to 
enhance drug loading, whilethe posterior side 
remained hydrophobic to better adhereto the 
capsular bag. Both lenses were drug loaded 
by soaking in solutions containing MXF and 
KTL, using conditions optimized in a previ- 
ous work [27] and sterilized by high hydrostat- 


ic pressure (HHP), a technique that proved 
to be adequate to sterilize hydrogels [28]. The 
IOLs were characterized with respect to some 
relevant physical properties and the adhesion 
of LECSs to the IOL surfaces was investigated. 
In vitro drug release tests were done and the 
antimicrobial activity of the released antibiot- 
ic was determined. Finally the in vivo perfor- 
mance of both lenses was compared. 


Experimental 


Disks with 1 mm of thickness and 5 mm 
of diameter and prototype IOLs of GF 
(Contamac, UK) and CI26Y (Physiol, Belgium) 
were used in the study. To increase the drug 
loading capacity of GF samples, their anterior 
surface was coated with CI26Y. Drug loading 
was done by soaking each diskIOL in a 
phosphate buffer solution (PBS) containing 
MXF and KTL (5 mg/mL of each drug). The 
disks/IOLs were sterilized in sealed bags 
containing the loading solution by HHP. The 
liquid uptake and the transmittance of the 
samples loaded with M XF+KTL and sterilized 
were assessed. 


In vitro drug release tests were performed 
in static sink conditions, by immersion of 
the disks in PBS. At specified times, ali- 
quots were collected and substituted by the 
same volume of fresh PBS solution, and the 
amount of drug released was quantified by 
UV-VIS spectrometry. 


The antibacterial activity of MXFKTL re- 
leased solutions, collected after 26 days of re- 
lease from sterilized disks, against S. aureus 
ATCC 25923 and S. epidermidis CECT 231, 
was investigated by agar diffusion tests. 


In vitro LEC adhesion tests were done by in- 
cubation for different periods: 24, 48 and 72 
h. Cells incubated into the wells without disks 
were used as control. After the different incu- 
bation times, the cells were stained with flu- 
orescein diacetate. The disks and the control 


were observed by fluorescence microscopy. 


In vivo experiments were done with Japanese 
white rabbits. Drug loaded and blank IOLs of 
both materials were implanted through cata- 
ract surgery in the right eye of the animals of 
the test and control groups, respectively. After 
the surgery, the usual eye drop therapy (Viga- 
mox 0.5 % and Acular 0.5 %) was applied to 
the right eyes of the control group. For thetest 
groups, no eye drops were instilled after the 
surgery. The eyes with implanted IOLs were 
inspected by slit-lamp microscopy, at pre-de- 
fined times between days land 21 Samples 
of aqueous humor were collected in the test 
groups on days 2, 7, 14 and 21and in the con- 
trol groups, only on days 2, 7 and 14, sinceno 
eye drops were administrated after the day 14. 


Results and discussion 


The CI26Y and GF samples loaded with M X- 
F+4KTL and sterilized by HHP presented lig- 
uid uptakes of 53% and 2%, respectively. The 
transmittance of both materials above 520 
nm assumed values >90 Y%. 


The in vitro drug release study shows that GF 
material release smaller amounts of M XF and 
KTL than CI126Y (Figure 1). In fact, GF disks 
present a lower liquid uptake, and therefore 
Shall incorporate a lower amount of drug. A 
controlled release of both drugs was achieved 
for more than 26 days for both materials. 
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Figure 1 MXF (A) and KTL (B) release profiles obtained from 
CI26Y and GF disks. 


(Adapted from J. Controlled Release 326 (2020) 245255) 


31 


32 


M icrobiological tests showed that the anti- 
bacterial activity of the released drug solu- 
tions against S. aureus and S. epidermidis 
was maintained after 26 days of release, as 
demonstrated by the presence of inhibi- 
tion halos (Figure 2). This confirms that the 
preparation and sterilization of the drug-load- 
ed materials did not induce activity loss of the 
antibiotic M XF. 


Figure2. Inhibition halos for drug solutions released atthe 26th 
day from Cl26Y (A) (10-fold dilution) and GF (B) Lfold dilution. 
Negative controls were obtained with sterile PBS. (Adapted 
from J. Controlled Release 326 (2020) 245-255) 


In Figure 3 it is possible to visualize an in- 
tense adhesion of LECSs onto the surface of 
both Cl26Y and GF disks, hydrated in PBS 
and sterilized by HHP, which increases with 
the incubation time. Incorporation of MXF 
and KTL reduced significantly the number of 
adherent cells, which became null for CI26Y. 


Figure 3. Fluorescent images of the LECs (green spots corre- 
spond to live cells) adhering to the surface of A) C126Y disks 
hydrated in PBS and loaded with MXF+ KTL and B) GF disks 
hydrated in PBS and loaded with M XF+KTL, after sterilization, 
for different incubation times. The control obtained without 
samples is also presented. (Taken from ). Controlled Release 
326 (2020) 245255) 


The performance of the studied IOL materials 
was finally assessed through in vivo tests. The 
concentrations of MXF and KTL in the aque- 
ous humour of the eyes of the rabbits of the 
test and the control groups are plotted as a 
function of time in Figure 4. MXF concentra- 
tions obtained with CI26Y lenses are always 
higher than for the correspondent control 
group, and stay above the minimum inhibito- 
ry concentrations (MICs) of S. aureus and 5. 
epidermidis for 21 days. In contrast, the con- 
centration of MXF measured in the control 
group is very close to the minimum values 
of the MICs. The concentration of KTL in the 
test group of CI26Y exceeds that found in the 
control group for the first seven days. In both 
groups (test and control), the concentration 
of KTL stays above the IC50 for COX-1 and 
COX-2 for 14 days that is the recommend- 
ed period for the instillation of NSAIDs eye 
drops after the cataracts' surgery. 


The results obtained with the coated GF IOLs 
are different due to the lower loading capac- 
ty of these devices, being the concentration 
of drugs released from these IOLs higher or 
similar to that measured in the control group, 
until day 7 for M XF, and for more than 2 days, 
but less than 7 days, for KTL. Even so, for the 
test group, the concentrations of M XF remain 
in the range of the MCs and of KTL above the 
IC50 of both COX for 21 days. 
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Figure 4. MXF (A and C) and KTL (B and D) concentrations in the aqueous humour after implantation of drug-loaded IOLs (º), and 
blank IOLs followed by eye drops therapy (º). Data for CI26Y are depicted in A and B and for GF in C and D. The minimum (dashed 
line) and maximum (full line) of the MICs for S. aureus (red) and S. epidermidis (yellow) (A), and of the IC50 for COX-1 (red) and 
COX-2 (yellow) (B) are shown. (Taken from J. Controlled Release 326 (2020) 245-255) 
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The eyes of the rabbits of the test and the con- 
trol groups were examined through slit-lamp 
examinations during the course of the ex 
periment (Figure 5). The mild inflammation 
observed in the test group on days land 3 
disappeared after day 7, meaning that the re- 


leased KTL was effective for the treatment of 
post-surgery inflammation. N either ocular in- 
fection/endophthalmitis nor PCO signs were 
detected in the eyes of the rabbits of both test 
and control groups. 


Days after the surgery 


Day 1 Day 3 


A - CIZ6Y 


Day 7 Day 14 Day 21 


ni eo] JD 


Blank IOLs 
+ Eye drops 


[616 


B - GF 


Coated 
drug-loaded 
IOLs 


Blank IOLs 
+ Eye drops 


eme! 
S 13673 


e 


Figure 5. Slitlamp photographs of the eyes from day 1 until day 21after surgery for: A) test groups with CI26Y IOLs drug-loaded (up) 
and control groups (down), and B) test groups with GF IOLs coated and drug-loaded (up) and control groups (down). (Taken from J. 


Controlled Release 326 (2020) 245-255). 


Conclusions 


Hydrophilic acrylic IOLs and hydrophobic 
acrylic IO Ls, coated with one hydrophilic layer, 
were successfully loaded with MXF and KTL 
to obtain controlled release platforms against 
ocular inflammation and endophthalmitis. 
Both types of devices seem to be more effec- 
tive than the commonly used eye drops thera- 
py because they ensure sustained concentra- 
tions of M XF and KTL in the aqueous humour, 


respectively above the MICs and the IC50 of 
both COXs, during the period recommended 
for the prevention of those ocular diseases. 
The drug-loaded IOLs kept adequate physical 
properties, and had no hazardous biological 
effects. Adhesion of LECs on the surface of 
the drug loaded IOLs was reduced in compar- 
ison with the blank IOLs, which may contrib- 
ute to reduce the PCO incidence. 
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ORTHOPEDIC BIOMECHANICS: COMPUTATIONAL 
ANALYSIS OF THE SHOULDER 


Carlos Quental (1), Madalena Antunes (1), 
João Folgado (1) 


1 /DMEC, Instituto Superior Técnico, Uni- 
versidade de Lisboa 


1. Introduction 


Biomechanics is the science that studies 
living organisms through the application of 
mechanical principles [1,2]. Dating its foun- 
dation back to the ancient Greece, with the 
series of Aristotle's essays De Motu Anima- 
lium (M ovement of Animals), in which ani- 
mals' bodies were studied as mechanical sys- 
tems to investigate the origin of movement, 
biomechanics owes its progress to many 
renowned scientists, including Leonardo da 
Vinci (1452-1519), Galileo (1564-1642), New- 
ton (1643-1727), and Borelli (1608-1769), who 
were instrumental for today's biomechanical 
concepts [3]. Contemporary biomechanics is 
a multidisciplinary field that is widely rec- 
ognized to play an important role in many 
different key areas [4]. Among its many ap- 
plication fields, we highlight the studying of 
the musculoskeletal system [5-7]; the cardio- 
vascular and respiratory systems [8-10]; the 
visual system [11,2]; and the multiscale in- 
fluence of forces on living organisms, which 
ranges from the macroscale (organs) to the 
nanoscale (cells) [13-15]. Regardless of its 
application, research in biomechanics has a 
main objective: to understand fundamental 
aspects of the human body. 


Over the past decades, the advance in com- 
puters and numerical methods promoted the 


computational modeling of biological sys- 
tems as a powerful and attractive solution for 
gaining insight into the human body biome- 
chanics and for addressing clinically relevant 
issues that are otherwise difficult, or even 
impossible, to investigate in-vitro or in-vivo 
[16,17]. Computational biomechanics can be 
broadly divided into two well-developed but 
separate modeling domains: multibody dy- 
namics and finite element modeling. M ulti- 


body dynamic models consider the complex 37 


interactions between human body segments, 
usually described as rigid bodies, to study hu- 
man movement [18]. Finite element models 
consider accurate representations of the hu- 
man body structures to account for localized 
tissue deformations [19]. 


Considering these two modeling domains, 
the aim ofthis paper is to showcase examples 
of our research in computational biomechan- 
ics, performed at the Institute of Mechani- 
cal Engineering (IDMEC), Instituto Superior 
Técnico, Universidade de Lisboa. For the sake 
of simplicity, all applications are related with 
the shoulder joint. For further applications, 
see our biomechanics group website: 

https:// biomec.idmec.tecnico.ulisboa.pt. 
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2. Applications in 
computational biomechanics 


a. Shoulder positioning durin 
superior capsular reconstruction 


The shoulder complex represents a per- 
fect balance between mobility and stability, 
which depends on the combined action of 
the rotator cuff— a group of 4 muscles, in- 
cluding the supraspinatus, subscapularis, 
infraspinatus, and teres minor— and other 
soft tissues. Due to the complex loading en- 
vironment that the rotator cuff is subjected 
to, rotator cuff tears (RCTs) are among the 
most common causes of shoulder disorders 
[20,21]. Recently, a novel technique, named 
superior capsular reconstruction (SCR), has 
been proposed to treat irreparable RCTs in 
younger, high-demand patients [22]. In this 
technique, the stability lost by an irrepara- 
ble RCT is compensated by the fixation of a 
graft between the superior rim of the glenoid 
and the humeral supraspinatus footprint. 
Although SCR has shown excellent short- 
term clinical outcomes, graft tear rates have 
ranged from 4.2% to 75% [23]. Considering 
that the position of the shoulder during the 
fixation of the graft may be a key factor for 
improving the outcomes of SCR, we have 
recently developed a multibody model of 
the shoulder to investigate the effect of the 
shoulder position during graft fixation on 
graft tear risk and shoulder stability after 
SCR for the treatment of an irreparable su- 
praspinatus tear [7]. 


The musculoskeletal model developed in- 
cludes 7 rigid bodies, representing the thorax, 
rib cage, clavicle, scapula, humerus, ulna, and 
radius, whose motions are constrained by 6 
joints; and are controlled by the main muscles 
of the upper limb. The SCR was simulated by 
deactivating the most common tor tendon 
of the rotator cuff (the supraspinatus tendon) 
and by adding a single-piece graft to the bio- 
mechanical model. The graft was modelled 
by 4 parallel segments, as illustrated in Figure 


1, and different shoulder positions of fixation 
were studied. The mechanical properties of 
the graft were defined based on experimental 
measurements performed at the Tissue Bio- 
mechanics Laboratory of Lisbon, located at 
IST [24]. 


Graft strain and shoulder joint reaction force, 
estimated through inverse dynamics [25], were 
used to evaluate graft integrity and shoulder 
stability respectively. Regarding graft integrity, 
shoulder positions during fixation above 20º 
of abduction presented a high risk of great 
tear when the arm was moved to its resting 
position, at the side of the trunk. For the re- 
maining shoulder positions of fixation (below 
20º of albduction), the results showed that 
SCR significantly improved shoulder stability 
compared with the preoperative condition; 
however, positions with the greatest improve- 
ments in stability were also the positions with 
the greatest risk for graft tear, due to their 
concerning larger graft strains. Motions in- 
volving excessive axial rotations for low arm 
forward flexions and for arm extensions be- 
hind the back were the most critical to the 
graft tear risk. Rehabilitation after SCR should 
be carefully designed considering these data 
to avoid compromising graft integrity Over- 
all, this study provided new and important 
information that assists orthopedic surgeons 
in the choice of the shoulder position for the 
fixation of the graft during SCR. 


ment 


egment 


Humerus 


Scapula 


(a) (b) 


Figure 1 Graft path for each segment in (a) anterior and (b) su- 
Es views, for a specific position of the shoulder during graft 
ixation. 


b. Biomechanical analysis of 
stemless shoulder implants 


For rheumatoid and osteoarthritic shoul- 
der joints, the shoulder arthroplasty, i.e, 
the replacement of the shoulder joint by an 
implant, has been shown to be a reliable 
treatment to reduce pain and restore shoul- 
der function [26]. Yet, it is associated with a 
set of complications that may compromise 
its long-term success [27]. The most con- 
cerning complication is aseptic loosening, 
whose etiology is multifactorial. Among the 
debated mechanisms of loosening, lack of 
primary stability and adverse bone adapta- 
tion due to stress shielding are key factors 
pointed for implant failure [28]. Primary 
stability is essential for osseointegration to 
occur, defined as the direct connection be- 
tween bone and an implant surface [29]. For 
a successful osseointegration, micromo- 
tions between bone and the implant are gen- 
erally accepted to have to be under 150 um 
[30]. According to Wollf's law, bone is a liv- 
ing tissue that adapts itself to the loads it is 
subjected to. When an implant is introduced 
into the bone, adverse changes in the me- 
chanical environment lead to negative bone 
adaptation processes that weaken the bone 
and, consequently, increase the risk of im- 
plant failure and periprosthetic fractures [31]. 


To avoid stem-related complications of tra- 
ditional stemmed implants, novel stem- 
less shoulder implants have recently been 
proposed. Although short- to mid-term 
outcomes show promising functional and 
radiological results, further investigation 
is necessary regarding their application as 
knowledge from stemmed shoulder im- 
plants cannot be transferred to these novel 
stemless implants due to their significant 
differences in design and rationale [28]. 
Considering five commercially available 
implants, we have recently studied the bio- 


mechanical performance of these implants 
regarding primary stability and bone adapta- 
tion using 3D finite element models [32,33]. 
The 3D geometry of a right humerus was 
reconstructed from the computer tomogra- 
phy data of the Visible Human Project [34] 
using ITK-SNAP [35]. Five stemless implants 
were modelled in Solidworks (Dassault Sys- 
têmes, Waltham, Massachusetts) based on 
the Eclipse Stemless Shoulder Prosthesis 
(Arthrex, Karlsfeld, Germany), Global Icon 
Stemless Shoulder System (DePuy, Warsaw, 
Indiana), SMR Stemless (LimaCorporate, 
San Daniele del Friuli, Italy), Simpliciti Shoul- 
der System (Tornier, Edina, Minnesota), and 
Sidus Stem-Free Shoulder System (Zimmer 
Biomet, Warsaw, Indiana). For the sake of 
simplicity these implants are referred to 
hereafter as “Eclipse-based,” “Global Icon- 
based,” “SMR-based,” “Simpliciti-based,” 
and “Sidus-based” models, respectively. 
The shoulder arthroplasties, depicted in Fig- 
ure 2, were virtually simulated in SolidWorks 
following the surgical guidelines available 
for each implant. 


(d) (e) 


Figure 2: Geometric models of the shoulder arthroplasties for 
the: (a) Eclipse-based, (b) Global Icon-based, (c) SMR-based, 
(d) Simpliciti-based, and (e) Sidus-based stemless implants. 
The green surfaces represent coated or blasted regions, while 
blue surfaces represent non-coated and non-blasted regions. 
The brown surfaces presented for the SM R-based, Simplici- 
ti-based, and Sidus-based implants represent the surfaces in- 
teracting with bone when contact was simulated between bone 
and the humeral head components. 
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The finite element models were developed in 
Abagus (Dassault Systêmes, Waltham, M as- 
sachusetts). Primary stability was studied for 
different bone densities through the evalua- 
tion of micromotions at the bone-implant 
interface for the osseointegration (coated or 
blasted) regions of the different stemless im- 
plants, and bone adaptation was evaluated 
by comparing the bone density distributions 
predicted for the intact bone model, with- 
out any implant, and the implanted models. 
Bone adaptation was simulated through the 
application of a mathematical model devel- 
oped within our research group, which mod- 
els bone adaptation as an optimization prob- 
lem that balances structural stiffness and a 
biological term related with subject-specific 
characteristics [36,37]. 


The Eclipse-based and Global Icon-based im- 
plants presented the best performance from 
the primary stability point of view, which may 
be due to their larger contact areas, com- 
pared with the remaining implants. However, 
regarding bone adaptation, the Global Icon- 
based design performed better, which sug- 
gests itto be a better overall solution among 
thesetwo designs. Figure 3 presents the abso- 


lute changes in bone density for each implant- 
ed bone with respect to the intact bone. Since 
stemless implants rely only on a metaphyseal 
fixation, their implantation is indicated only 
under good bone quality conditions, for which 
all stemless designs presented micromotions 
below 150 um. Considering this, other stem- 
less implants might also be better long-term 
solutions. From the bone adaptation point 
of view, the SM R-based model presented the 
best performance. Another interesting finding 
of our studies was that the contact between 
the humeral head implant and the resected 
bone surface, recommended by some man- 
ufacturers, had 2 opposite outcomes: on the 
one hand, it provided greater primary stabili- 
ty but on the other hand, it negatively impact- 
ed bone adaptation. 


Although further investigation is necessary 
to prove the benefits of stemless implants, 
they are currently being considered the fu- 
ture of shoulder arthroplasty [38]. Computa- 
tional biomechanics is a powerful tool not 
only to predict their long-term performance, 
but also to improve implant design and im- 
plantation technique. 
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Figure 3: Absolute changes in bone density Au between the implanted models and the intact bone model for the fully completed con- 
dition of osseointegration for the: (a) Eclipse-based, (b) Global Icon-based, (c) and (f) SMR-based, (d) and (9) Simpliciti-based, and 
(e) and (h) Sidus-based stemless implants. The results in (c), (d), and (e) were obtained considering contact between the humeral 
heads of the implants and the resected surface of the bone, whereas the results in (f), (9), and (h) were obtained considering no 
contact. The anatomical terms lateral, posterior, medial, and anterior are denoted by L, P, M, and A, respectively Bone resorption is 
depicted in red (A <-0,2 g.cm3) and orange (-0.2 g.cm'3 <Au <-0.19.cm'3), whereas bone apposition is depicted in light blue (0.1 
g.cm3>Au >0.2 g.cm3) and dark blue (Ay >0.2 g.cm3). White denotes an equilibrium condition (-0.1g.cm'3 <Ay <0.1g.cm'3). 


3. Concluding remarks 


Computational modeling of biological sys- 
tems has been shown to be a valuable tool 
to examine both physiological functions and 
clinical problems. In this paper, research in 
orthopedic biomechanics performed at ID- 
MEC was showcased to demonstrate how 
computational biomechanics can assist clin- 
ical decision-making and contribute to im- 
proving the overall quality of medical treat- 
ments and of people's well-being. 


Despite the significant advances in compu- 
tational biomechanics over the last decades, 
many challenges remain. Considering its ul- 
timate aim of translating results to the real 
world (e.g., the clinical setting), two main 
challenges must be addressed. Validation of 
computational biomechanical models is es- 
sential to convince the medical and scientif- 
ic communities of their utility. This requires 
confirming that biomechanical models are 
a trustworthy representation of a biological 
system relative to precise measurements 
of analogous real-world systems, which is 
extremely challenging due to the complexi- 
ty and variability of biological systems. The 
other challenge regards the development of 
subject-specific models. Most biomechani- 


cal studies in the literature are based on a 
limited number of anatomical datasets, and 
results are extrapolated to the general pop- 
ulation. For a wider applicability in the clin- 
ical setting, the inclusion of subject-specific 
characteristics in computational models is 
crucial. Moreover, because the development 
of subject-specific biomechanical models in- 
volves technically- and time-demanding pro- 
cedures, novel accurate, time-effective, and 
cost-effective methods are mandatory. Sta- 
tistical shape modeling and machine learn- 
ing frameworks have recently demonstrated 
high potential for tackling this challenge. 


In Portugal, biomechanics has made con- 
siderable progress in recent years, demon- 
strated by the international recognition of its 
different research groups working in biome- 
chanics. We have the expertise to be at the 
forefront in this area and to contribute to the 
solution of challenges that will have signifi- 
cant societal impact. 
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BLACK HOLES AND RIEMANN-HILBERT PROBLEMS 


M. Cristina Câmara and Gabriel Lopes Cardoso 
Center for Mathematical Analysis, Geometry and Dynamical Systems, 
Department of M athematics, Instituto Superior Técnico, 
Universidade de Lisboa, Av. Rovisco Pais, 
1049-001 Lisboa, Portugal 


1. Introduction 


The Einstein field equations of General Rel- 
ativity are a system of 10 coupled second or- 
der non-linear partial differential equations 
(PDE's) in 4 independent variables for the 
space-time metric g,w for a given distribution 
of matter and energy in space-time. 


Despite their elegant form, these PDE's are, 
in general, very complicated to solve. How- 
ever, exact solutions to these PDE's may be 
found under simplifying assumptions, for 
instance by imposing symmetries, such as 
spherical symmetry. Exact solutions to the 
Einstein field equations play an important 
role, since they model important gravita- 
tional phenomena, such as black holes. One 
exact solution that is of particular math- 
ematical and physical importance is the 
Schwarzschild metric. It describes a static 
black hole in space-time, and it constitutes 
the first exact solution to the Einstein field 
equations that was obtained [1]. 


It turns out that rather than directly solving 
the above system of non-linear PDE's, ex 
act solutions to the Einstein field equations 
in vacuum (and, more generally to the field 
equations of some gravitational theories 
without a cosmological constant) can be ob- 
tained by means of a Riemann-Hilbert (RH) 
approach to the gravitational field equations. 
In this approach one obtains exact solutions 
by solving a matricial Rigmann-Hilbert prob- 
lem (RHP). 


What then is a RHP? RH problems can be 
formulated in simple terms, as follows. Con- 
sider first a function & that is analytic in two 
regions D*and D- ofthe complex plane which 
are separated by a common boundary, such 
as the real line or the unit circle. Then, one 
may determine the jump of & across T by 
studying the boundary values of & that one 
obtains by approaching the boundary T from 
each side. 


A RHP is the inverse problem: given a jump 
across T, one wants to determine a function 
that is analytic in both regions D* and D 
and that satisfies the given jump condition. 
Therefore, a RHP is the problem of recon- 
structing an analytic function on CAT from 
its jump across T. À simple example there- 
of is as follows: suppose that one wants to 
determine two functions & and &, analyt- 
ic and bounded in D* and D, respectively, 
and continuous on F suchthat & =& on 
F. Then the unique solution is d, = D=c, 
where C denotes a constant, as a conse- 
quence of Liouville's theorem in complex 
analysis. 


It is a remarkable fact that a big variety of 
problems in mathematics, applied mathe- 
matics and physics can be rephrased as a 
RHP. In doing so, there are always two steps 
involved. The first step consists in reducing 
the original problem to a RHP. This means 
that the original problem, formulated in a 
very different setting, is reduced to a problem 
in complex analysis. This often means, as in 
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the case of the gravitational field equations, 
replacing the original non-linear problem by a 
linear one. The second step consists in solv- 
ing the RHP. Here, one should note that the 
level of difficulty in solving scalar valued and 
matrix valued RH problems differs greatly: 
the latter ones are much harder to solve. The 
RHP associated to solving the gravitational 
field equations is a matricial RHP. 


2. Lax pair 


As mentioned above, exact solutions to the 
gravitational field equations can often be 
found by restricting to solutions that possess 
symmetries. In the following, we will consider 
the subspace of solutions that possess suffi- 
clently many commuting isometries so that 
the gravitational field equations become re- 
duced to a system of coupled second order 
non-linear PDE's in two spacetime dimen- 
sions, depending on two coordinates, (p,v), 
called Weyl coordinates, which take values in 
the Weyl upper half plane, i.e. pe R* and v E 
R. In two dimensions, the PDE's to be solved 
take the form 


d(pxA)=0, (1) 


where A is a matrix one-form, 


A=M-tdM M=M(p,v), (2) 
with components 4 = 4, do +A, dv. Here M 
= Mé is a matrix that is invariant under an 
involution *, often referred to as generalized 
transposition. This involution is determined 
by the gravitational theory that one is study- 
ing. In (1) * denotes the H odge star operator 
in two dimensions, which is defined by 


x*dp = —o dv 


*dvo = dp. (3) 


where o = +], depending on whether the Weyl 
upper half plane is endowed with an Euclide- 
an (o=1) or with a Minkowski metric (o=- 1). 
In the case of the Einstein field equations in 
vacuum, the involution % is matrix transposi- 
tion, and the matrixM is a two-by-two matrix, 


ira tá + B2/A da) | (4) 


B/A A 


Here, A and B are functions of the Weyl coor- 
dinates, and together with another function 
Wy, they determine the space-time metric Buy 
The function y is obtained from M by inte- 
gration, 


Oy = 


zP Tr (AS = o A2) , 
Lo Tr(AÇA,) (5) 


Ow = 


Therefore, if one solves the non-linear PDE's 
(1), one will have obtained a solution g,, of 
the gravitational field equations. The question 
then is how to determine a matrix M(p,v) that 
solves (1). This can be done as follows. 


The non-linear PDE's (1) form an integrable 
system. Following [2] a system of (non-lin- 
ear) PDE's is an integrable system if it can 
be represented as a compatibility condition 
of an auxiliary overdetermined linear system 
of differential equations, called a Lax pair 
for the given non-linear system. A complex 
parameter must be included in the Lax pair. 
In the case of (1), the Lax pair is the Breiten- 
lohner-M aison linear system for the unknown 
matrix X [3], involving a complex parameter 7, 
called the spectral parameter, 


T(dX + AX) =+dX. (6) 


The spectral parameter 7 is not a constant: 
It must satisfy an algebraic relation involving 
another complex variable w as well as the 
Weyl coordinates p and v. This relation, called 
the spectral curve, is defined by 


By solving this relation for 7 as a function of 
w, p and v, one obtains the following two 
functions, 


vI+rop (8) 


Now it is important to note the following: let 
us pick a function qo ofthe form (8) and con- 
sider the associated Lax pair (6) with a given 
one-form 4 as an input. Then the Lax pair (6) 
Is a linear PDE for the unknown X. However, 
since our goal is to solve (1) in order to obtain 
a solution M, from which we then infer the 
one-form 4 = Mº dM, one must actually look 
at the Lax pair (6) as an equation for the pair 
of unknowns (X,4). Therefore, at first sight, it 
appear that one has not gained anything by 
introducing the Lax pair (6). 


But this is actually not so. By making use of 
a specific RH problem that will be introduced 
momentarily one simultaneously obtains a 
solution M(p,v) to the field equations (1) as 
well as a solution X (T=q, (pv), pv) to the 
Lax equation (6) with input 4=M" dM, pro- 
vided the RHP can be solved explicitly. The 
RHP is question consists in obtaining a ca- 
nonical bounded Wiener-H opf factorization 
of a monodromy matrix, with respect to an 
appropriate contour in the complex plane 1. 
Let us introduce this particular RHP. 


3 Canonical factorization of a 
monodromy matrix 


Let M (cv) be a matrix function of the complex 
variable «w and denote by M the matrix 
that is obtained from M(«w by composition 
using the spectral curve AR (7), 1.e. 


ST), 6) 


Mol T) = M (: T 


defined foral TEC (0), pe R*, ve Rsuch 
that M,, (D) Is invertible. The matrix (9) Is 
referred to as monodromy matrix 


In (9), we consider (p,v) as an arbitrary pair 
of parameters in a neighbourhood of (p,,v,), 
and take 7 as the (complex) independent vari- 
able. To emphasize this aspect, we have de- 
noted the matrix Min (9) by Ms (T). 


To proceed, we need to make a few assump- 
tions [4]. The most important ones are as 
follows: 


Assumption 1: There exists an open set S in the 
Weyl upper-half plane such that, for every (p.v) € 
S, one can find a simple closed curve T in the 7- 


plane, satisfying 


1. T passes through the fixed points of the trans- 


formation T -—s —£:; 

É ; 

2. T encircles 7 = 0; 
3. àfrTreT then -SeET; 


Mp, 
to-invariant open set O in the 
taining TV, and such that Mult 
on O. 


v(T) and its inverse are analytic in an 
T-plane con- 


)=Mpu(T) 


4] 
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Let D* denote the simply connected interior 
regionof= and D=C(D*'U T). 


Assumption 2: MpulT) admits a canonical 
bounded Wiener-Hopf factorization with respect to 
P, that is, a matriz factorization of the form 
Mou(T) = Mp v(T)X(7, p, v), (10) 
where X(T.p.v) (respectively MovlT )) and its 
inverse are analytic and bounded in the 
set D* UO (respectively D”" U O). 
X(T.p.v) satisfies the 
X(O, p,v) = Inxn. 


open 
Moreover, 


normalization condition 


Note that if M,, (7) admits such a canonical 
factorization, then it is unique. 


Next, let us define 


(11) 


pu 


M(p,v) = Jim Mor). 


Under the above assumptions, one proves 
the following theorem [4]; 


Theorem 1 M(p,v) is a solution of the field 
and X(T = 


Polp.v).p.v) is a solution to the Laz pair (6). 


equations (1) on the open set S, 


The proof of this remarkable theorem is based 
on the solution of a second Riemann-Hilbert 
problem and highlights the deep role of the 
spectral curve, the normalization condition in 
the factorization and the choice of the con- 
tour [4]. 


É 


Figure 1: The exterior region of the Schwarzschild black hole 
solution [5]. 


4. An example 


Let us consider a solution of the Einstein field 
equations of General Relativity namely the 
Schwarzschild black hole, and show how it 
can be obtained in a simple manner by per- 
forming the canonical factorization of the fol- 
lowing two-by-two matrix, 


U— th () 


0 q Sir 


w—m 


M(w) = | - (12) 


where the parameter m is taken to be posi- 
tive. We consider the case o=1 (the case o=- 
1 can be treated analogously). 

Substituting «w by the expression on the right 
hand side of the spectral curve relation (7), 
we obtain 


(T—Ta)(T+l/T9) 


(r—Ti)(r+1/71) 


(13) 


where, 


[ , 
v-m>—y(m—v) + p? 


T(Pp.v) = 
b 


[ , 
vem—y(m+ vv) + pº 


p (14) 


T2lp.v) 


Now let us pick a contour T'such that r, and 7, 
do notlieon 7. There arefour possible classes 
of contours that are determined by the posi- 
tions of 7, and 7, with respect to the contour. 
Let us pick a contour such that 7, and 7, are 
inside the contour 7. Performing the canoni- 
cal factorization of (13) with respect to 7” we 
obtain, 


=m-(T)ma(7T) (15) 


with 


E) 
+ 
“ 


m(T) = 


= 
o tw 
s Es 


= 


E Pe 
q 


m-AT) = 


(16) 


Note that m, is analytic in D* and satisfies 
the normalization condition m, (0)=1, while 
m is analytic and bounded in D: Then, using 
(11), we obtain 


Õ mo) (00) 


(17) 


This is of the form (4) with À =7, /, and 
B=0. 
The resulting spacetime metric g, , reads 


Adte + AH 
+A pêdg?, 


“ (do? + dv”) 
(18) 


ds4 = 


with e” obtained by straightforward integra- 
tion of (5). Finally, using the bijection from 
the domain r> 2m, 0 < 60 <r in spherical 
coordinates onto the Weyl upper half-plane, 
given by 


p = vrº-2mr sind, 


vV = (W- 


(19) 


m) cost, 


A becomes 4 = 1- 2m/r, and the space-time 
metric g,, takes the form 


A f > ) 
ds = (1 E Jae + (1 +] dr? 
r ro 


r2 (dO? + sin? 0 dó”) . 


(20) 


which is the familiar form of metric describing 
the exterior region of the Schwarzschild black 
hole (see Figure 1), obtained here by means 
of the canonical factorization (15) using sim- 
ple complex analytic tools. 


Finally, note that by choosing a contour be- 
longing to any of the other classes of con- 
tours and proceeding in a similar manner, we 
would obtain other distinct space-time solu- 
tions to the Einstein field equations [4]. 
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In the last decade we have seen a truly rev- 
olutionary development in the study of grav- 
itation. Gravitational waves offer an entirely 
new channel of discovery, and allow to un- 
derstand two-body systems as never before. 
The GRAVITY instrument - where Técni- 
co plays an important role - is now able to 
monitor flares close to the supermassive 
black hole at the center of the M ilky Way. The 
Event Horizon Telescope, a virtual, Earth 
sized device, has produced the first imag- 
es of such object. There is much more to 
come in the next decades. The experimen- 
tal breakthroughs rely on, but also encour- 
age and motivate theoretical work in strong 
field gravity. With the hope of giving an im- 
pression of what we do, here we present a 
collage of the work being performed in the 
Técnico gravity group grit /CENTRA. 


|. Introduction 


General Relativity (GR) is a monumentto the 
human intellect, and our best description of 
the gravitational interaction. N evertheless, it 
seems to contain its own demise: it predicts 
the existence of black holes (BHs), in the 
interior of which the classical theory fails, 
and for which we still don't have a quantum 
completion. We need to collect data to un- 
derstand where and how GR fails. Do black 
holes even exist, and are they the objects 
predicted by GR? 


Luckily GR is a geometric theory of gravity 
which admits physical wavelike solutions, a 
feature which is imposed by a finite maximum 
speed of propagation in our universe. These 
gravitational waves (GW) areto gravity what 
light is to electromagnetism. Their physical 
effect is to stretch and compress spacetime 
transverse to their direction of propagation. 
Unfortunately for experimental work, they are 
extremely weak. For the purposes of detec- 
tion, sufficiently strong GWs can be generat- 
ed only by pairs of extremely compact astro- 
physical objects. Examples of such systems 
are pairs of BHs or neutron stars. In fact, the 
first direct detection of GWs was made onlyin 
205 by LIGO, one full century after Einstein 
produced the theory itself. The first detection 
turned out'to be emitted from the collision of 
two BHs with around thirty times the mass of 
the Sun [1]. By the very nature of BHs, which 
do not emit light, this system could only be 
seen through GW'S. In the seven years since 
then, many more signals have been detected. 
The detected waves all agree perfectly with the 
predictions of GR. Thus, although the search 
to find an experimental discrepancy goes on, 
for now GR remains the flagship description 
of gravity We must therefore take the lessons 
of the theory seriously. 


Obtaining and understanding solutions to 
the gravitational field equations is thus an im- 
portant problem in its own right. In contrast 


51 


52 


to those of Newtonian gravity however, the 
field equations of GR are nonlinear. Despite 
their geometric beauty, they are moreover, 
highly complicated, involving many unknown 
variables that interact in a subtle manner. Con- 
sequently, finding exact solutions in generic 
physical situations is a hopeless task. Various 
approaches to obtain approximate solutions 
have been developed. These include pertur- 
bative setups, such as post-N ewtonian, linear 
and nonlinear perturbation theory [2, 3), but 
also the use of rigorous mathematical esti- 
mates, pure computational approaches [4-6] 
and indeed modeling. Phenomenological 
studies also shed light on mysterious aspects 
of the solution space. In this short article we 
give an overview of the research being per- 
formed inthe grit /CENTRA (gravitation in 
Técnico) group [7]. 


Il. Extreme spacetimes 


In highly dynamical regions of spacetime, 
where we can appeal to no small parameter 
for linearization, approximate analytic meth- 
ods necessarily fail to give accurate predic- 
tions, and rigorous mathematical estimates 
of solutions are typically too abstract in na- 
ture to draw physical conclusions. In such 
situations we turn to numerical relativity: the 
construction of approximate solutions of Ein- 
stein's equations by computational means. 
The compact binary spacetimes discussed 
in the introduction constitute an important 
example, but from a purely theoretical point 
of view there are, perhaps, even more funda- 
mental possibilíties. 


A foundational issue concerns the failures of 
GR: we know BHs cloak infinities (so-called 
spacetime singularities), is it possible that 
some form “naked” to our eye, in such a way 
that we can see and study them? One can 
think, loosely, of a naked singularity as being 
a region in which the spacetime geometry it- 
self “blows-up” in a way that observers watch- 
Ing from arbitrarily far away can feel the blow- 
up, rather than it being hidden inside a BH. 


The weak cosmic censorship conjecture says 
that given (non-special) initial data for space- 
time, naked singularities should not arise as 
we evolve to the future. Mathematical tools 
are far from establishing the conjecture in 
generality but it can be tested numerically by 
evolving members of families of initial data 
that interpolate between “weak” and “strong. 
The idea is that if we could find (open) sets 
of initial data on which naked singularities do 
form, then we would have a convincing coun- 
terexample to the conjecture. 


Numerical experiments of this type have been 
performed since the 1990s [8, 9], albeit with 
the majority assuming spherical symmetry. 
There are numerous subtleties and details, 
but the striking result is that, in spherical 
symmetry there is compelling evidence that 
(for typical matter models) there is a unique 
threshold solution that lies between disper- 
sion and collapse to a BH. In other words, 
spacetime either eventually settles down and 
becomes peaceful, or collapses to form a BH, 
with the resulting singularity safely hidden 
from far-away observers. Precise tuning to the 
threshold of BH formation results in a unique 
spacetime with a naked singularity and with a 
remarkable, simple, mathematical structure. 


The open question is to what extent this struc- 
ture is maintained when the solutions have 
less symmetry In this setting we may ask 
what happens to GWSs as they collapse under 
their own self-gravity Numerical evolutions 
without symmetry are more computationally 
intensive, and so the answer has remained 
a mystery for thirty years as both numerical 
methods and availability of large scale com- 
pute clusters has improved. Nevertheless, a 
growing body of work [10, 11, 13-16], to which 
members of grit have made a major contri- 
bution, indicates that although aspects of the 
spherical phenomenology survive, the gen- 
eral story is much more complicated. Thus, 
despite this progress, a complete picture of 
the threshold of collapse, for now, remains 
elusive. In Fig. 1we plot the gravitational field 


(a quantity called Kretschmann curvature sca- 
lar) ffom one of our spacetime simulations 
close to the threshold of collapse. There one 
can see fine features forming in the curvature. 


FIG. 1 A snapshot of the log of the (Kretschmann) curvature 
scalar in a vacuum spacetime close to threshold of BH forma- 
tion (see [10, 11] for more details). The spacetime was computed 
numerically using the bamps numerical dia code [2] de- 
veloped, within a larger collaboration, at grit.T 


each separately whilst added communicating data be- 


is axially symmetric. The large peak visible lies on the axis of 
symmetry. 


II. Waves and the global 
structure of spacetime 


Astrophysical sources lie millions of light- 
years away. Thus, between our detectors and 
the source, there are trillions of wavelengths 
of a GW. Conceptually it's like having the de- 
tectors infinitely far away from the source. 
How do we capture such a large distance 
numerically? The solution came through con- 
formal compactification, a mathematical con- 
cept brought to GR by Roger Penrose (The 
2020 Nobel in Physics, see [17). Conformal 
compactifications are not new; in fact sailors 
and cartographers have used them for centu- 
ries to map the surface of the Earth onto a 
piece of paper: stereographical projection. In 
this transformation, the north/south-pole of 
the terrestrial sphere maps to infinity in the 
plane. This transformation changes areas but 
keeps angles the same (shapes) — that is why 


Antarctica looks so large in worldmaps! By 
means of a conformal transformation of the 
spacetime metric in GR, we can bring “in- 
finity” to a finite location in the transformed 
spacetime. However, in GR there is more 
than one type of infinity as depicted in the 
Penrose-Carter diagrams in Fig. 2. There is 
spatial infinity iº representing infinity reached 
along spacelike geodesic curves —such 
as the solid red lines in a); future and past 
null-infinity j*, denoting infinity along null di- 
rections — straight lines in b) depict the path 
of light; and future and past time-like infinity 
i£, representing infinity along time-like geode- 
sic curves — like the dashed orange lines in 
a). These notions are not just a mathematical 
curiosity. In fact, several physical quantities 
such as the mass and angular momentum 
are only rigorously defined at infinity. More- 
over, fundamental questions of the theory 
such as the cosmic censorship conjecture 
can be phrased in terms of properties of the 
global structure of spacetimes. 


Conformal compactification is a method that 
can be applied to the Einstein equations to 
reach any infinity [18]. However, of special rel- 
evance is future null-infinity 7*, because it is 
there that GWs are unambiguously defined. 
The complex astrophysical scenarios that cre- 
ate GWSs requires the use of numerical relativ- 
ity which in tum relies on slicing spacetime 
to solve the equations. Among the possible 
spacetime slicings, there are two promising 
options that reach future null-infinity One is 
CCE (Cauchy-characteristic evolution [19]) — 
see c) in Fig. 2, where inner standard space- 
like slices (in red) are matched to outer null 
slices (in blue). An alternative is the use of 
hyperboloidal foliations [20]: space-like slices 
that do not end at iº, but rather rise up and 
intersect null-nfinity as shown in purple in c). 
The end goal of these efforts is to numerically 
evolve and read the gravitational waveforms 
at j* directly This requires the synergy of in- 
vestigations in mathematical relativity as well 
as the development of robust numerical tech- 
niques which are currently being pursued in 
QE 
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a) 


r= const 
— t= const 


FIG. 2. Penrose-Carter diagrams for the Minkowski spacetime. 
The different elements of the conformal boundary are shown, 
along with different slicings probing these regions. 


IV. Gravitational wave 
phenomenology 


We discussed how to extract GWSs after they 
propagate through spacetime, but what pro- 
duces them? A typical source is a binary of 
two BHs orbiting each other, moving in an in- 
ward spiral (the inspiral) as they lose energy 
due to the emission of GWSs, until they col- 
lide and form a larger BH. After this merger 
stage, the final BH relaxes to a stable state, in 
a process that is called ringdown. During it, 
the BH vibrates according to some character- 
istic frequencies named quasinormal modes 
[21]. This three-stage behaviour is very well 
understood and using a combination of ana- 
Iytical techniques and the power of numerical 
relativity we are able to predict how the GW 
signal of two BHs looks like. This knowledge 
is crucial because not only does it contribute 
to the detection of the GWSs, but also allows 
us to understand the properties of source 
populations, such as the mass and angular 
momentum of the BHs, their number densi- 
ty distance to us, etc. 


u=const 
v=const 


GWSs are produced in the most extreme re- 
gions of spacetime, and are an unparalleled 
tool to test GR. By comparing the theoretical- 
ly predicted waveforms with the ones mea- 
sured, we can test GR. 


But GWSs also have the potential to probe new 
physics that we know exists. Around 30% of 
the Universe is composed by matter that does 
not emit light - dark matter. There is a vast 
zoo of options for its fundamental nature, 
each one leading to new GW features. For 
example, some models suggest dark matter 
is composed by new ultralight particles with 
masses of — 102º eV (compared to these, an 
electron has a gigantic mass of 0.5 x 10º ev). 
If they exist, they can grow around spinning 
BHs and condensate into clouds through a 
phenomena called superradiance [22]. These 
clouds, which have a structure that resembles 
that of the hydrogen atom, emit monochro- 
matic GWs which are already being searched 
by current detectors [23] Quite poetically, 
we are doing particle physics with classical 


gravity! Dark matter and other types of exotic 
matter can form ultracompact objects almost 
as compact as BHs but that do not possess 
an horizon from which light cannot escape. 
These objects have a rich phenomenology 
with perhaps their most striking feature being 
the presence of echoes in the ringdown stage. 
They correspond to repetitions of the original 
signal that is being reflected back and forth in 
the gravitational potential [24], see Fig. 3 By 
searching them, we can infer if we are actual- 
ly observing BHs or BH-impostors. GWs can 
also tell us about the environment where they 
were created. Binaries evolving in galactic 
centers are surrounded by an accretion disk 
and dark matter that generate torques. The 
latter can drag the binary, accelerating its in- 
spiral. Different accretion models then result 
in different waveforms [25, 26]. Finally, we can 
also detect GWs coming from far-away dis- 
tances that were originated in the very early 
Universe and may contain precious informa- 
tion on its formation and evolution. Atgrit, 
we are lucky to have several world leading ex- 
perts working in some of the problems listed 
above. 


GW strain 


0 50 100 150 
time [ms] 


FIG. 3. The “signal” when stars falls onto a large BH (dashed 
black line) or onto an exotic compact object, a bBH impostor 
(solid red). Echoes are a distinctive feature of new physics at 
the horizon. 


V. The future of 
gravitational physics 


Now is the right moment to study BHs and 
regions with violent spacetime dynamics: the 
future of gravitational physics promises to 
be extremely exciting. LIGO and Virgo, now 
joined by the Japanese detector KAGRA, are 
expected to restart operations in early 2023 
with increased sensitivities and are planned 
to reach their design sensitivity in a couple of 
years [27]. On the longer-term, a next gener- 
ation of ground-based observatories, includ- 
ing the Einstein Telescope [28] in Europe and 
Cosmic Explorer in the USA [29], are being 
planned [30]. These detectors will be signifi- 
cantly more sensitive than LIGO and are ex 
pected to come online in the 2030's. Ground- 
based detectors are especially sensitive to 
GWs emitted by BH binaries with tens of 
solar masses. Detecting signals from much 
more massive BHs, which emit at much low- 
er frequencies, requires going into space. 
This is what LISA [31], a space-based GW ob- 
servatory led by the European Space Agency, 
promises to do. LISA has been approved to be 
launched in the next decade and will observe 
GWSs emitted by supermassive BH s with mil- 
lions of solar masses or by extreme mass-ra- 
tio-inspirals, systems in which a small com- 
pact object with tens of solar masses orbits 
around a supermassive BH, such as the one 
in the center of our own galaxy. At even lower 
frequencies, pulsar timing arrays 
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FIG. 4. Characteristic strain versus frequency for a variety of de- 
tectors and sources. The area between the detector's curve and 
a given a source indicates how loud that signal is. Created using 
gwplotter [42] See Ref. [43] for details. 
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(PTASs) are expected to detect in the next few 
years the ensemble of GWs emitted by bina- 
ries of supermassive BHs [32]. The sensitivity 
band of these different detectors and the pro- 
totypical sources they can detect are illustrat- 
ed in Fig. 4. All these observatories will allow 
us to increase the precision of current mea- 
surements but also increase the sheer num- 
ber of GW detections. This will open an era 
of precision and data-driven GW science, that 
promises to start a new era of discoveries in 
the gravitational physics. 


The observation of GWs has also given us 
the opportunity to make the first joint de- 
tection of the gravitational and electromag- 
netic waves emitted by a binary neutron star 
merger, which gave birth to a new form of 
multi-messenger astronomy [33]. The future 
of this new form of astronomy promises to be 
bright and will allow us to learn more about 
the astrophysics of compact objects [34], but 
also test whether GWs propagate differently 
than electromagnetic waves [35-37] and pro- 
vide independent measurements of how fast 
our Universe is expanding [38] 


While the community was witnessing the rise 
of GW astronomy, the Event Horizon Tele- 
scope collaboration captured the first direct 
images of a supermassive BH [39] including 
the recently revealed image of SgraA * [40] the 
supermassive BH at the centre of our own 
galaxy. Yet, this is justthe beginning of BH im- 
aging. The community is now putting in place 
the science case for a next generation Event 
Horizon Telescope (ngEHT) with increased 
sensitivity [41]. 


Given these developments on the observa- 
tional side, in the grit group we are actively 
preparing for the promised discoveries that lie 
ahead. We are members of the LISA consor- 
tium, of the Einstein Telescope collaboration 
and are leading a working group on funda- 
mental physics for the ngEHT collaboration. 


VI. Conclusion 


In this short article we presented a broad pic- 
ture of the research we are carrying out in the 
grit group. Our main goal is simple: we 
want to understand gravitation and useit as a 
tool to learn more about our Universe. Given 
that observations of GWs and BHs are final- 
ly allowing us to compare observational data 
with the predictions of GR in the presence of 
the most extreme gravitational fields in the 
Universe, this has become even more press- 
ing. Achieving this goal is no easy task, but we 
hope to have convinced you that it is certainly 
worth pursuing. 
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